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We construct initial data for several black holes with arbitrary momenta and spins by 
a new method that is based on a compactification of Brill-Lindquist wormholes. When 
treated numerically, the method leads to a significant simplification over the conventional 
approach which is based on throats and an isometry condition. 



In this paper we wilLuse the standard formulation of the ADM initial value 
problem as given by YorkEJ. 

We will use g a b to represent the physical 3-metric, which will be related by a 
conformal factor ip to the flat 3-metric by the relation g a b = ip^gab- Likewise, the 
physical extrinsic curvature of the spacetime is denoted by K a b and is related to the 
conformal curvature by K a b = ip~ 2 K a b. Furthermore, we use the relation K a a = 0. 

1 Brill and Lindquist Data 

The Brill and Lindquist initial data model is given by 

There are certain advantages and disadvantages to using this in a numerical code. 
Advantages: 

1. There is no need for an elliptic solve and no possibility of numerical error. 

2. No need to worry about existence or uniqueness of the solution. 
Disadvantages: 

1. There is a coordinate singularity at f = ri n \. Numerically, however, this presents 
surprisingly little difficulty. Data with this sort of singularity present has been 
evolved previously in 3D codea3'u. Furthermore, it is expected that numerical codes 
of the future will use apparent horizon boundary conditions to remove regions of 
the grid inside apparent horizons, and that would naturally exclude these points. 

2. The very thing which made this problem so simple, namely time-symmetry 
Kat = 0, severely restricts the set of spacetimes we can construct. Specifically, we 
can neither impart boosts nor spins to the holes. 



2 The Method of Throats and Images 

The strategy for this procedure is to excise spherical regions from the conformally 
flat grid, and use an isometry condition to supply the boundary conditions on the 
spherical surfaces. 
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The extrinsic curvature is given by a base form, which has boost and spin 
parameters. One adds together a number of these base forms equal to the number 
of black holes in the initial value problem to obtain an extrinsic curvature with the 
appropriate ADM values on each hole (though these parameters can only truly be 
identified with spin and boost in the limit of wide separation). 

The base solution to the momentum constraint for a single black hole is: 

T^ab / nQ„& i Tjb^a ( ~ab ^ab\ nc„ \ , / acd o ~ ~~b , bed o ~ 

K = ^-^ yF n + F n — [g — n n ) F n c ) H — - (e b c ridn + e b c ridn ) 

~'~ ' (2) 

where r is the radius from the black hole's center, n a is the normalized radius vector 
for r, P a is the boost vector for the hole, and S a is the spin vector. This base form 
of the extrinsic curvature does not yet satisfy the isometry condition at the throats, 
so we must now add an infinite number of well-defined "image terms" together to 
accomplish this. This method has several advantages, and has been implemented 
in a variety of settings by Cook et. alia As before, we compare advantages and 
disadvantages of this method. 
Advantages: 

1. No coordinate singularities - the spherical regions which contain them have been 
excised from the grid. 

2. Kab ^ 0, and thus boosts and spins can be given to each black hole separately 
(in some sense) 

Disadvantages: 

1. The method is difficult to code, inner boundary conditions present a complication 
which break most "off-the-shelf" elliptic solvers. 

2. It is more difficult to construct numerical approximations when infinite series 
are involved. It is difficult to say anything analytically about these data sets. 

3. Whether solutions exist and are unique is generally unknown (except for a few 
casesH) . 



3 A New Method 

We propose a new methoofl which, essentially, extends the Brill and Lindquist 
solution by using the base momentum constraint solution used in the method of 
throats and images. The solution is of this form: 



~, 1 = E 9 |^ | > P=Uk ab k«» (3) 
— ' 2 \r — r(„) | 8 



(n) 



With this ansatz the Hamiltonian constraint can be written as an elliptic equation 
for u on all of R 3 (no singularities at the r^), and one can use a generic elliptic 
solver to obtain a solution for u, which is regular over the entire grid. 
Advantages: 

1. Can accommodate boosts and spins. 

2. Simple to code. 

3. Less difficult to approximate. We can construct solutions that are first order 
accurate in (3 - the boost/spin term - for single black holes by expanding u as 
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u = l + eui + ... and writing e(3 instead of just (3 in Eq.(|3|). Note that this first order 
correction is an upper bound to the full non-linear solution. One can easily use this 
to check whether a given numerical solution makes sense. Any correct numerical 
solution must have the property 1 < u < u\ everywhere on the grid. In one of 
the plots below, a first order approximate analytic solution is compared to a fully 
non-linear numerically generated solution for a single boosted and spinning black 
hole. 

4. Existence and uniqueness can be provenQ. We can also determine that u has no 
minima. Further, if comparing two solutions that differ only in the overall magni- 
tude of the boost/spin term (3, the solution for u with the larger boost/spin will be 
everywhere greater then the solution with the smaller boost/spin. 
Disadvantages: 

1. While there no longer are singularities at the punctures, u is only twice diffcr- 
cntiable there. With the standard second order stencils, convergence is not second 
order near the punctures, which however is confined to a region close to the punc- 
tures. The tiny dip at the punctures in figs. 1 and 2 is a genuine feature, as can be 
seen for higher resolution in both the numerics and the approximation. 



Analytic Approximations 

Initial data for two black holes 
A2BS data set on z— axis 




Figure 1: Convergence study for spin- Figure 2: Comparison between analytic 
ning black holes undergoing spiral infall approximation and full numerical solve 

for a single boosted spinning black hole 
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